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Abstract 

We show that any weak solution to the full Navier-Stokes-Fourier system ema- 
nating from the data belonging to the Sobolev space IT 3 ' 2 remains regular as long as 
the velocity gradient is bounded. The proof is based on the weak-strong uniqueness 
property and parabolic a priori estimates for the local strong solutions. 
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1 Introduction 

In continuum mechanics, the motion of a general compressible, viscous, and heat con- 
ducting fluid is described by the thermostatic variables - the mass density g = g(t, x) and 
the absolute temperature $(t, x), and the velocity field u = u(t, x) evaluated at the time 
t and the spatial position x belonging to the reference physical domain Q C R 3 . The time 
evolution of the fluid, emanating from the initial data 

q(0, ■) = go, 0(0, •) = O , u(0, ■) = u in Q, (1.1) 

is governed by the Navier-Stokes-Fourier system of partial differential equations: 
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Equation of continuity: 



d t g + diy^u) = 0; 
Momentum equation: 



;i.2) 



d t (gu) + div x (gu g> u) + V z p(^>, #) = div x S($, V x u); 
Total energy balance: 



;i.3) 



-£>|u| 2 + £>e(£, i?) ) u 



div x (p(^,^)u) (1.4) 



div x (s(tf, V x u)u) - div x q(tf, V x $). 



The symbols p = pig,^) and e = e(g,$) stand for the pressure and the (specific) 
internal energy, respectively. Furthermore, S = S($, V x u) denotes the viscous stress given 
by 

Newton's rheological law: 



(tf, Vsii) = //(tf ) V^u + V^u - -div.ui , 



;i.5) 



while q = q($, V^tf) is the heat flux determined by 



Fourier's law: 



q(l?,V 8 u) = -K(iJ)V I i>. 



;i.6) 
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We restrict ourselves to bounded regular domains Q C R 3 , with energetically insulated 
boundary. More specifically, we consider the standard 



NO SLIP BOUNDARY CONDITION: 

u\an = 0, (1.7) 

together with 




where n is the (outer) normal vector to dQ. Under these circumstance, the total mass 
as well as the total energy of the system are constants of motion, specifically, integrating 
equations (11. 2p . (jl.4p over Q we obtain 

g(t, ■) dx= [ Q dx = M , (1.9) 
n Jn 

f n (^l u ! 2 + £ e (£' ^)) (*» ') dx = J Q (^o|u | 2 + Q e{g , t? )) dx = E . (1.10) 

From the mathematical viewpoint, the Navier-Stokes- Fourier system suffers the same 
deficiency as most of its counterparts in continuum mechanics - the lack of sufficiently 
strong a priori bounds. It is known that the problem (11.11 - II. 8p is well posed, locally 
in time, in the framework of classical solutions, see Tani |26j, or in slightly more general 
energy spaces of Sobolev type, see Valli [27], [28J, Valli and Zajaczkowski, [291], among 
others. Besides these local existence results, Matsumura and Nishida established in a series 
of papers [20], [21], [22] global- in-time existence of strong/classical solutions provided the 
initial data are sufficiently close to an equilibrium solution. Similarly to the well known 
incompressible Navier-Stokes system, the existence of global-in-time smooth solutions for 
any physically admissible and large initial data remains an outstanding open problem. 

The system (11.11 - IT78T) has been also studied in the framework of weak solutions. Hoff 
and Jenssen [T4] established global existence for radially symmetric data in R 3 . They 
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also identified one of the main stumbling blocks in the analysis of the Navier-Stokes- 
Fourier system, namely the (hypothetical) appearance of vacuum zones, where the density 
vanishes and the classical understanding of the equations breaks down. More recently, 
Bresch and Desjardins [I], [5] discovered a new a priori bound on the density gradient 
leading to global-in-time existence in the truly 3D— setting conditioned, unfortunately, 
by a very specific relation satisfied by the density dependent viscosity coefficients and a 
rather unrealistic formula for the pressure that has to be infinite (negative) for g — > 0. 

1.1 Weak and dissipative solutions 

We adopt the approach to the Navier-Stokes-Fourier system originated in [9] and further 
developed and detailed in [TT]. Suppose that p and e are interrelated through 

GlBBS' EQUATION: 

>&Ds{g, •&) = De(g, tf) + p(g, $)D , (1.11) 

where s = s(g,$) is a new thermodynamic function called (specific) entropy. If g, are 
smooth and bounded below away from zero and if u is smooth, then the total energy 
balance ( 11.41) can be replaced by 

Thermal energy balance: 



Furthermore, dividing (11.121) on d we arrive at 

Entropy (production) equation: 



with 



d t {gs(g, #)) + div x (gs(g, $)u) + div 2 



'q(fl,Vsfl)' 
i9 



(7, 



Entropy production rate: 



;i.i3) 



a 



d 



S(i?, V x u) : V x .u - 



1.14) 



Let us remark that the systems §U\ Ol OD , (Ol Ol H32D , and ([H Ol [LTD are 
perfectly equivalent in the class of classical solutions. 
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Another crucial observation, which is the platform of the approach developed in |llj . 
is that we can relax (11.141) to 

o > I (s(0, V.u) : V.u - SL y^l (1-15) 

provided the system is supplemented by the integrated total energy balance ( 11.1011 , mean- 
ing the relations (11.21 11.31 II. 10111.131 . with ( 11.151) are still equivalent to the original system 
(11.21 11.31 H-41) . at least for smooth solutions. The resulting problem is mathematically 
tractable, more specifically, we report the following results: 

• Existence. The problem { (Ol-IO). (Ol- fTTOD . (031) . (051 } admits global- 
in-time weak solutions for any finite energy initial data under certain structural 
restrictions imposed on p, e, s, and /x, k presented in Section [2] below, see [TTl 
Theorem 3.1]. 

• Compatibility. Any weak solution of { (Ol-IO). (TL51 - ITTDjl . IQ3| . } 
that is regular solves the original problem (11.11 - 11.80 , see [Hi Chapter 2] . 

• Weak-strong uniqueness. Any weak solution coincides with a classical solution 
emanating from the same initial data provided the latter exists, see [T2]. 

1.1.1 Relative entropy and dissipative solutions 

The property of weak-strong uniqueness mentioned above is closely related to the relative 
entropy inequality introduced in [12]. We start by introducing another thermodynamic 
function: 

Ballistic free energy: 
H e ( e ,&) = ge{g,#) - egs(g^), (1.16) 

see [7], together with 

Relative entropy functional: 



i ff |u - U| 2 + H e (g, t?) - dH& ^ 9) (g - r) - H e (r, 6) 



dx, 
(1.17) 

where {g, i9, u} is a weak solution of the Navier- Stokes- Fourier system and {r, G, u} is an 
arbitrary trio of smooth functions satisfying 



r > 0, > 0, U 



<->n 



0. 



'1.181 
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As observed in [12j, the weak solutions of the Navier-Stokes- Fourier system satisfy 

Relative entropy inequality: 



E ( g , A u|r, O, U)] Z + /; 1 1 (m V..) = V.. - ■»<* V f - V ' tf ) dxdt (1.19) 

- X ~ U ^ ' 5 * U + ~ u ) ® u : V * U - P(6> ^) div * u ) dx dt 

+ f f (s(0, V x u) : V^U dx dt 
Jo Jn v 

-ft (q(s(q,$) - s(r,Q))d t Q + g(s(g,$) - s(r,Q))u-V x Q) dx dt 

J « aZ 

JO Jf7 77 

+ Jo L ( ( X ~ r ) ^ P(r ' 6) ~ r U ' Va;P(r ' 9) ) dX dt 
for any trio {r, G, U} satisfying (j!.18p . Taking {r, 0, U} the strong solution of the same 
system emanating from the same initial data gives rise to a Gronwall type inequality for 
£ yielding g = r, # = 6, u = U, see [T2] . 

As the proof of the weak-strong uniqueness principle uses only the integral inequality 
(I1.19P without any reference to the original system of equations, we may introduce a new 
class of dissipative solutions satisfying solely the relative entropy inequality (11.191) for any 
trio of admissible test functions {r, G, U}. Note that a similar concept was introduced by 
DiPerna and Lions [19] in the context of the Euler system. 



1.2 Regularity criteria 

A conditional regularity criterion is a condition which, if satisfied by a weak solution, 
implies that the latter is regular. Similarly, such a condition may be applied to guarantee 
that a local (strong) solution can be extended to a given time interval. The most celebrated 
conditional regularity criteria are due to Prodi and Serrin [23], [21] and Beal, Kato and 
Majda [3], Constantin and Fefferman f6] for solutions to the incompressible Navier-Stokes 
and Euler systems. Recently, similar conditions were obtained also in the context of 
compressible barotropic fluids and the full Navier-Stokes-Fourier system, see Fan, Jiang 
and Ou [8], Sun et al. [25], and the references cited therein. 

In view of the results of Hoff et al. [13], [15], certain discontinuities imposed through 
the initial data in the compressible Navier-Stokes system propagate in time. In other 
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words, unlike its incompressible counterpart, the hyperbolic-parabolic compressible Navier- 
Stokes system does not enjoy the smoothing property typical for purely parabolic equa- 
tions. Analogously, a solution of the full Navier- Stokes- Fourier system can be regular only 
if regularity is enforced by a proper choice of the initial data. 

Our approach to conditional regularity is based on the concept of weak (dissipative) 
solutions satisfying the relative entropy inequality ( II. 19ft : 

• In Section [21 we introduce the structural restrictions imposed on the thermodynamic 
functions p, e, s and the transport coefficients /i, k so that the Navier- Stokes- Fourier 
system may possess a global in time weak (dissipative) solution for any finite energy 
initial data. 

• Now, suppose that the initial data are more regular, specifically, 

o < q < go, o < < O , go, O e w 3 ' 2 (n), u g w 3 > 2 (yi- r 3 ), (1.20) 

satisfying, in addition, the necessary compatibility conditions 

u |an = 0, V I O ■ n| an = 0, V a .p(^ o ,0o)|an = div rE §(0 o ,u o )|an. (1.21) 

Under these circumstance, the Navier-Stokes-Fouries system admits a local-in-time 
strong solution {g, 0, u} constructed by Valli [28] such that 

g, G C{[0,f};W 3 ' 2 {Q)), G L 2 (0, f; W 4 ' 2 (Q)) , d t $ G L 2 (0, T; W 2,2 (Q)), (1.22) 

u G C([0, f ]■ W 3 ' 2 {tt- R 3 )) n L 2 (0, f; W 4 ' 2 (n; R 3 )), d t u G L 2 (0, f; W 2 ' 2 (n ; R 3 )). 

(1.23) 

• Modifying slightly the arguments of [12J to accommodate the strong solutions be- 
longing to the regularity class (11.22)1 we show in Section[3]that the associated (global- 
in-time) weak solution {g, 0, u} emanating from the initial data coincides with the 
strong solution {g, 0, u} on its existence interval [0,T]. 

• We show that if 

limsup (\\g(t, -)|U 3 . 2 (fi) + ||^(*, -)\\w^(n) + ||u(t, ^Ww^iQ;^)) = oo, (1.24) 
t->f- 

then, necessarily, 

limsup ||V x u(t, -)IU°°(n;Ji 3x3 ) = oo, (1-25) 
t-tf- 

see Section HI 
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• Finally, in Section we observe that the solutions constructed by Valli are classical 
(all necessary derivatives are continuous) for any t > 0. We conclude by the follow- 
ing conditional regularity result for the weak solutions of the Navier- Stokes- Fourier 
system: 

Suppose that {g, u} is a weak (dissipative) solution of the Navier-Stokes-Fourier 
system in (0, T) x Q, emanating from the initial data belonging to the regularity class 
specified through (11.201) . (jl.2ip . and satisfying 



Then u is a classical solution of the Navier-Stokes-Fourier system, unique in the class 
of weak (dissipative) solutions, 

see Theorem 12.11 below. 

2 Preliminaries, main results 

We start with a list of structural restrictions imposed on the functions p, e, s, and 
k. The interested reader may consult [l~Tj Chapter 1] for the physical background and 
possible relaxations. We suppose that the pressure is given by the formula 



ess sup || V x u(t, ■) 11x00(0.^3x3) < oo. 
te(o,T) 




(2.1) 



with P G C^O, oo) n C 3 (0, oo) satisfying 



P(0) = 0, P'(Z) > 0, < 



\P{Z) 



P'{Z)Z 



< c for all Z > 



(2.2) 



Z 




(2.3) 



Accordingly, in agreement with Gibbs' equation (11. lip . 




(2.4) 
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and 

a ^) = s{^)+~ (2-5) 

where 

gW = -| fP(Z) / (Z)Z ,fc^=0. (2.6) 

Finally, the transport coefficients are continuously different iable functions of the tem- 
perature satisfying 

Ml +^ A ) < M#) < A*(l +^ A ), < c for a11 ^ e [0,oo), - < A < 1, (2.7) 



k(1 + i? 3 ) < k(0) < 7S(1 + tf 3 ) for all 7? G [0, 00). (2.8) 
2.1 Weak and dissipative solutions 

It was shown in (TTJ Theorem 3.1] that under the hypotheses (12.11 - I2.6I) . the problem { 
( ll.ll - ll.3p . ( ll.lOp . (11.131 - IT7T5T) } admits a global- in-time weak solution for any initial data 
satisfying 

g , ^0 e go, i? > a.a. in Q, u G L 2 (fi; i? 3 ). (2.9) 

Moreover, the weak solution {g,$,u} satisfies the relative entropy inequality (11.191) . see 

rjj]. 

Accordingly, a trio {g,$, u} will be called dissipative solution of the Navier-Stokes- 
Fourier system (11.11 - 11.81) provided it obeys the relative entropy ( 11.191) for any choice 
of smooth functions {r, 0, U} satisfying ( 11.181) . 

Summarizing the results of [HI Chapter 3] and [12] we obtain: 

Proposition 2.1 Let Q C R 3 be a bounded domain of class C 2+u . Suppose that the 
thermodynamic functions p, e, s and the transport coefficients fi, n obey the structural 
hypotheses A2.1\ - {2Jty . Finally, let the initial data belong to the class specified in A2.9\) . 

Then the Navier- Stokes- Fourier system possesses a dissipative solution {g, $, u} on an 
arbitrary time interval (0, T) that enjoys the following regularity properties: 

g>0a.a. zn(0,T)xfi, g G C([0, T]; L 1 ^)) n L°°(0, T; L 5/3 (fi)) n lf((0, T) x Q) (2.10) 
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for a certain (3 > |; 

tf>0 a.a. m (0, T) x SI, d G L°°(0, T; L 4 (fi)) n L 2 (0, T; V^ 1,2 (fi)), (2.11) 

^bg^el^O.r;^)); (2.12) 

ueL 2 (0,T;W«(n;R 3 )), « = ^ ^ C weak (0, T; L 5 / 4 (fi; ,R 3 )). (2.13) 

2.2 Main result 

In accordance with the programme delineated in the introductory part, our main goal is 
to show that for regular data the weak solution remains regular as long as we can control 
the gradient of the velocity field. More specifically, we will show the following result. 



Theorem 2.1 Under the hypotheses of Proposition l2~l\ let {g, t?, u} be a dissipative 
solution of the Navier-Stokes-Fourier system on the time interval (0, T) belonging 
to the regularity class 112. 1Q{ - 12.13\) , with the (regular) initial data satisfying U.20\) . 
together with the compatibility conditions $1.21]) . 

Suppose, in addition, that 

ess sup HVajU^, ■) 11x00(0.^3x3) < oo. (2.14) 
*e(o,T) 



Then {g, $, u} is a classical solution of the Navier-Stokes-Fourier system satisfying 
(TO-rOj) in (0,T) x tt. 



Remark: For technical reasons, we have omitted the effect of the bulk viscosity in 
the viscous stress S(?9, V x u). The possibility to extend the result to more general forms of 
the viscous stress is discussed in Section^ 

Here, classical solution means that all functions and all derivatives appearing in the 
equations (11.21 - I1.4D are continuous in (0, T) x Q, the functions g, u, together with 
their first order derivatives, are continuous in [0, T] x Q and satisfy the initial conditions 
(11.11) as well as the boundary conditions (11.71) . (11.81) . 
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Although the rest of the paper is essentially devoted to the proof of Theorem 12.11 we 
will obtain a few other results that may be of independent interest. 



3 Local existence and weak-strong uniqueness 

To begin, it is convenient to introduce a scaled temperature 



5 = KOd), where K0&) = / k(z) dz. 

Jo 

Accordingly, the thermal energy balance fll . 12j) reads 

q de(g,{}) 



K{$) 



S(0, V x u) : V.u 



u • V X E) - A: 

$ dp(g,i}) 
K(d) dd 



(3-1) 



div^u 



3.1 First a priori bounds 

Our goal is to derive a lower and upper bounds on the temperature. We first rewrite (13. ip 
in the form 

(d t E + u • V X E) - DAE = DAE - Bdiv x u E, 

where 



D 



i -i 



,4 



S(i?, V^u) : V x u, B = — t — s — — — ( q- 



K{$) 



K(d) d$ 



Now, in view of the hypotheses (12. ip . (12.41) . and (12.81) . it is easy to check that the exist 
constants B_, B such that 

< B < B(t, x) < B for all t, x. 
Applying the standard comparison argument, we therefore deduce that 



S(r,-) > inf H(0,x)exp (-B / ||div x Ti(t, 0IU«(n) dt ), r > 0. (3.2) 

xGO. \ Jo J 

In order to obtain an upper bound, we need a similar estimate for the density, namely 

(3.3) 

12 



inf g (x)exp - / Hdiv^u^, -)IU~(n) dt 

x£U \ Jo 



L°°(Q) 



(If 



< q(t, •) < sup Qo(x) exp ( f \\div x u(t, 

that follows easily from the equation of continuity (II. 2ft . 
Now, we may use hypotheses (12.71 - 12.81) to observe that 

A(t r ) <A\\V x u(t,-) lll-^xs); 

whence 

S(r,-)< (3.4) 



sup 5(0, x 



exp (d A ||V x u(t,-)||ioo(n.ji3xs) dtj + exp (b J ((div^i, ■) \\ L °°(n) dt\ 
for r G [0, T], where D depends only on 

fT 

inf go(x), inf S(0, x), and / div x ii £«>(m di. 
x£fl xefi Jo 

3.2 Local existence of strong solutions 

Rewriting the system of equations (11.21) . ( II. 3p . (11.121) in terms of the new unknowns 
{g, S, u} and taking the a priori bounds (13.2)) . (13. 4p into account, we can apply the local 
existence result of Valli [281 Theorem A and Remark 3.3]. Going back to the original 
variables, we obtain: 

Proposition 3.1 Under the hypotheses of Proposition \2J\ suppose that the initial data 
{£ 0) ^cb u o} satisfy U.20\) . and the compatibility conditions M.21\) . 

Then there exists a positive time T such that the Navier-Stokes-Fourier system U.l\ - 
\1.8\) admits a unique strong solution {g, u} in the class 

g,de C([0, f}; W 3 ^ 2 (n)), u G C([0, f ]; W^{^ R 3 )), (3.5) 

& G L 2 (0, T; W 4,2 (£l)), d t $e L 2 (0,f;W 2 > 2 ({l)), (3.6) 
u G L 2 (0, T; W A ' 2 (Vt] R 3 )), d t ue L 2 (0,f;W 2 ' 2 (n;R 3 )). (3.7) 
Moreover, there exist scalar functions g, g, i? ; depending solely on 

rf 

inf g Q , sup^o, inf^oj supf? , and on / || V z u|||oo m .R3x3) dt, 
n n n n Jo v . ; 

suc/i t/iat 

< g(r) < q(t, •) < g(r), < g(r) < 0(r, •) < #(r) (3.8) 

for any r G [0, T]. 
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3.3 Weak strong uniqueness 



We claim that in view of the result [12] and its generalizations obtained in [10J, the dissipa- 
tive solution obtained in Proposition 12.11 coincides with the strong solution of Proposition 
13. li on the time interval [0, T] provided they start from the same initial data (11.201) . (11.211) . 
This deserves some comments since both [T2] and [TU] deal with classical solutions having 
all relevant derivatives continuous and bounded in (0,T) x Q. 

Given the integrability properties of the dissipative solutions stated in (12.101 - 12. 1 3[) 
and the regularity of the strong solutions (13.51 - I3.7[) . we can easily check that the trio 
{r = g, 6 = U = u} can be taken as test functions in the relative entropy inequality 

dm. 

Now, following step by step the arguments of [101 Section 6, formula (79)] we deduce 
from (I1.19P the estimate 



S (^g, u g, u 



t=T 



t=0 



< 



+ c 







I _ ~ II 2 _L 
l U U llwl>«(n;il3) + 



d - 



dt 



(3.9) 



+ 



X2( 



Xi(t)£ (g,&,u 
14 



g, u) dt 



u — u 



1 + g + q\s(q, #)| 

where, similarly to [ID] , we have denoted 

h = [h] css + [h] TCS = h- [h] ess , 

[h} ess = *{Q,&)h, $ G C c °°((0,oo) 2 ),0 < $ < 1, 
$ = 1 in an open neighborhood of a compact K C (0, oo) 2 
where K is chosen to contain the range of [g, specifically, 

[£(t, x) , &(t, x)} G # for all x G H, t G [0, T] . 

The functions Xj, i = 1,2 are of the form 

Xi(t) 

= bi(t)(l + ||Qt^||ioo((o,f)xn) 4" II V 2 $|| T.aoffn fwni + H9tu|| r.oorcn fwm + l|V 2 u| 



| dxdt, 



lL°°((0,T)xQ) 



lL°°((0,T)xfi) 



l°°((o,t)xq) y J 



where &, are bounded positive functions determined in terms of the amplitude of g, u 
and their spatial gradients in [0, T] x Q. Focusing on the most difficult term, we have 



Xi(t) / e\ s (e,$)\\ l u -u] rcs | dx < xz(t) \\[qs(q 



, v)\res ||L4/3(Q) 



u 



u 



L 4 (f2;_R3) 



14 



^ II ~ i i 2 

< e u — u 



< e ||u — u||wl,a(Q._ R 3) res 1 1 £4/3 (q) 

g, d, u) for any e > 0, 



i w i, a(Q;i? 3) + c(»x20)£ 3/2 (e, u 

where, exactly as in [TDl Section 6.1], we have used the structural properties of s and the 
embedding 

W l ' a (ty M> L 4 , a = E < A < l. 



5-A' 5 



Going back to (13. 9p we may infer that 



£ [g, $,11 g, u 



+ c 
*=o ./O 



I _ ~ II 2 _l_ 

l U U llw 1 '° ! (n;i? 3 ) + 



< 



0-0 

0, u) d£, 



2 



(3.10) 



where, by virtue of (13. 61) . (13. 7ft . 

Applying Gronwall's lemma we obtain the following conclusion: 



Proposition 3.2 Under the hypotheses of Proposition lSAi let {g, 0, u} be the local strong 
solution specified in Proposition ^. 1\ and {g, 0, u} a dissipative solution obtained in Propo- 
sition \2.1\ emanating from the same initial data. 
Then 

g — g, = 0, u = u in [0, f] x fi. 



4 Conditional regularity 

Our goal is to show that the energy norm 

\\g(t, -)||h/ 3 . 2 (q) + ||0(*, -)llw 3 ' 2 (n) + ||u(t, •) 11^3,2(^3) 

associated to a strong solution of the Navier-Stokes-Fourier system remains bounded in 
[0, T] as long as 

SUp ||Va;U(t, ■)IU° (fi;B 3x3 ) < G < OO. (4.1) 

te(o,T) 

We remark that we already know that 



< q(t) < g(r, ■) < g(r), < < 0(r, •) < 



(4.2) 



for any r G [0,T], see ( 13. 8ft . where the bounds depend only on G, the initial data, and 
the length of the time interval. 
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4.1 Energy bounds, temperature 

Multiplying equation (13.11) on 



q <9e(£,$) 
we obtain, after a routine manipulation, 



AS, 



ess sup \\E(t, -)\\ w i,2 {n) + 
te(o,T) 



L 2 (0,T;L 2 (n)) 



+ 



I < — 1 1 1 2 



\L 2 (0,T;W 2 - 2 (n)) 



<c(S,data). (4.3) 



Since § is already known to be bounded, the estimates (14.31) transfer to Indeed note 



that 



Did 



-i 



K~\E) Dl~ + 



k-He) 



where, since H satisfies the homogeneous Neumann boundary condition, the term V X S 
may be estimates in terms of the Gagliardo-Nirenberg inequality 



|VxS|| L 4 (n . fl 3) < iihUloo^iiasii^^). 



(4.4) 



Consequently, relation (14.31) implies 



ess sup \\${t,-)\\ w ia {n) + 

te(o,f) 



L 2 (0,T;L 2 (Q)) 



+ 



t) 



L 2 (0,T;W 2 > 2 (O.)) 



<c(5,data). (4.5) 



4.2 Energy bounds, velocity 

Our next goal is to deduce similar bounds for the velocity field. To this end, we write the 
momentum equation in the form 

a,u - WMl V.6) = -u ■ V.fi - i^Vj + IMMlv.J (4.6) 

£ Q OV Q OQ 

ldp{g,&) _ 
= hi + - — V x £, 

where, in accordance with the previous estimates 

sup || hi (t,-) \W(n-,RS) < c(B, data). 
te[o,f] 
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Taking the scalar product of (14. 6 j) with — div^t?, V x u) we obtain 

S(tf, u) : d t V x u dx + f \ div x .§(^, V^u) 2 dx 
<> Jn g 



where, furthermore, 



hi + \ dp ^ d \ x g \ ■ div x §($, V^u) dx, 
n v q OQ I 



S(#, u) : d t V x u dx 



d y fifi) 



dt Jn 4 



V x u + V'.u — -div-ul 



a 4 



V x ii + Vlu div x ul 



dx 



dx. 



In view of the estimate (14. 5 j) we conclude that 



d y /x(^) 



dt 4 



V x u + V;u div^ul 



dx + 



div x S(tf, V^u) dx (4.7) 



where 



hi + i ^f^ Vsg) ■ div x §(tf, V x u) dx + f h 2 dx, 
n \ g oq I Jn 



h 2\\ L ^o,f;L\n) < c(S,data). 



4.3 Elliptic estimates 

In order to exploit (14.71) . we have to show that the L 2 -norm of div :r S(^, V x u) is "equiva- 
lent" to the L 2 -norm of V 2 u. We have 



H = div x §(tf, V^u) = n(&) ^Au + -V^div^fi 
where, in accordance with (14.51) . 

l|h3|lL°°(o,T;L 2 (Q;i? 3 ) — c(B, data). 



+ h 



3; 



Consequently, we may integrate (14. 7p with respect to t and use the standard elliptic 
estimates to obtain 



ML^n-flS) dt < c(£,data) 1 + / \\V x q\\ ) > T e [°> T ]- ( 4 - 
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Now, we differentiate the equation of continuity with respect to x to obtain 

dt (d Xi g) + u ■ V x (d Xi g) = -d Xi (u) ■ V x g - {d Xi g) div x u - ^.div^u; (4.9) 

whence 

— / \V x g\ 2 dx < c(B, data) [ (\V x g\ 2 + \V x g\\V x div x u\) dx. (4.10) 
dt Jq Jq v ' 

Combining (14.81) . (14.101) with the standard Gronwall type argument we get the follow- 
ing estimates: 

sup || V x g(t, •) \\l 2 (q-,r3) < c(B, data). (4.11) 

t£[0,f] 

Next, due to (jOjl . 

ll^^llL 2 (0,f;L 2 (n ;j R3)) + ll^llL 2 (0,T;H/ 2 > 2 (n;H 3 )) — C{B, d&t&), (4-12) 

and, finally, by means of the equation of continuity, 

II^Hl 2 (o,T;L 2 (c)) < c(B,data). (4.13) 

4.4 Energy estimates for the time derivatives 

Differentiating the momentum equation (11.31) with respect to t and setting V = d t u we 
obtain 



g (d t V + u • V.VJ - div x S(#, V,V) = - (d t gV + gV • V.uJ + hi, (4.14) 

with 



2 



n'ifydtd) ( v x u + v*u - ~di V:c ui 



-V, (d t p(gj)), 



hi = -d t gu ■ V x u + div x . 
where, in accordance with the previous estimates, 

ll h ilL 2 (o,f;Ty-i.2 ( n;-R3)) < c(G, data). (4.15) 

Seeing that 



- (dtgV + gV ■ V,u) ■ V = \V\ 2 gdiv x ii 



we may take the scalar product of (I4.14p with V and, integrating over Q, we deduce the 
energy estimates for V = d t u: 

sup ||9tu||L2(n ; fl3) + ||^tu||ii2(o,f.wi,a(n;flS)) < c(G,data). (4.16) 

t€[0,T] 
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and, according to Sobolev's embedding theorem, 

11^11^(0,^6(0;^)) < c(G, data). (4.17) 
Going back to (11. 3p we compute 

(Au+iv.div^u) (4.18) 

= m* + u • V.u) - (v.u + V .V - ^div.ul) V J + VMC, *); 
which, combined with (14.161) and the previous estimates, gives rise to 

sup \\u(t, •) \\ W 2, 2(q ;R 3) < c(G, data). (4.19) 
te[o,f] 

Next, bootstraping (14.181) via the elliptic regularity, yields 

o ||u(<, 11^(0;^) dt<c(l + J Q T \\V x g(t, Wl^rs) dtj . (4.20) 
Furthermore, multiplying (14. 9 j) on | V^^V^ yields 

^I|V x fi»||x,8 ( n ; jia) < c (ll^'^lll 6 (o ;j R3) + ||u||wa.«(n ; jis)||V a .e||i6(n.jis)) • (4-21) 
Thus, combining (14.201) . (I4.2ip . we may infer that 

sup [HV^t, OIU 6 ^ 3 ) + \\dtQ(t,-)\\LS(n;R3)} + \\u\\ L 2(of-w^(n;R3)) ^ c(G, data). (4.22) 

te[o,f] 

4.5 L p — estimates for the temperature 

Our next goal is to apply the technique of L p — L q estimates to the parabolic equation 
(I1.12p . For this purpose, we first need Holder continuity for i?. To this end, let us write 
equation fll . 12j) in terms of e(g, E) = e(g, K _1 (S)). 

gd t e(g, E) + gu ■ V x e(g, E) - A X E = V x u) : V x u - p(g, tf)div x u := h. 

Note that 

A X E = div x ( V » e , ) - div, ( d ^ E) V x g) . 

\e e (g,E)J \d e e(g,E) J 
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We obtain 



V x e \ , , J; _ ( d g e { g,E) 



gd t e(g, E) + gu- V x e(g, E) - div x - — ~ = h + div x - ' ~ V x g , 

\d e e{g,^)J \d@e(g,^) J 

which, after dividing on both sides by g, yields 

d t e(g,E) + (u- — ^|-~r ) • V x e(g,E) - div x ^' l< 1 



g 2 d e e(g,E)J ~ \gd e e(g,E)J 

= -+ dA ~^l IV^I 2 + div. ( dA ~ 6 ^ V x g] . (4.23) 
g g 2 d e e(g,E) \gd e e(g,E) J 

Note that V ' x ~g G L°°(0, T; L 6 (Q, R 3 )) according to ffl~22l) . we can apply the standard 
theory of parabolic equations with bounded measurable coefficients to deduce that 

e(g, E) is Holder continuous in [0,T] x Q. 

See Ladyzhenskaya et al. [TB] . Since p is already Holder continuous in the set [0, T] x f2 
(cf. the estimates (14.221) ). we find 

E(hence$) is Holder continuous in [0, T] x £7. (4.24) 

Now we write (If .121) in the following form: 

d t E + u ■ V X E - DA X E = Dh. (4.25) 



with the diffusion coefficient D = ( g de( ^f^ ) being Holder continuous. The LP — L q 



-l 

>1\ 

theory for parabolic equations is now applicable yielding 

G L p (0, f; W 2 ' 6 (Vl)), dtti G L p (0, f; L 6 (fi)) for any 1< p < oo. (4.26) 

see Amann [1] , [2] , Krylov [16] . Note that the integrability in the spatial variable is limited 
by the integrability of the initial data. 

Remark: Applying similar arguments to the momentum equation (11.31) . we could 
obtain analogous estimates for the velocity field: 

u G L p (0,f; W 2 '\VL; R 3 )), d t u G L p (0,f; L 6 (fi; R 3 )) for any 1< p < oo, (4.27) 

however, we do not need this refinement in the future analysis. 
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4.6 Full regularity 

Differentiating (14.251) with respect to time yields 

d t E t - DA x E t = (Dh) t + D t A x E - (u ■ V x E) t : = h. (4.28) 
According to fl4TTB|) . fT4~22l and (g2ZD, it is easy to check that 

ll^ll^(o,f^(n)) <c(G,data). (4.29) 
Standard energy method and elliptic estimates yield 

ll^llL°°(o,f;L 2 (n)) + II— illi^o.iw^Cfi)) — c(G , data) , (4.30) 

ll-llL°°(o,7w 2 . 2 (n)) + ll"llL 2 (o,f ; iy 3 . 2 (n)) ^ c(G, data), (4-31) 

as well as 

II— tllL^CCTW 1 . 2 ^)) + ll ,=; «llL 2 (o,f';L 2 (n)) — c(G,data), (4.32) 

II— llL°°(o,f;W 3 ' 2 (n)) + ll-<llL 2 (o,f ; H' 2 ' 2 (n)) — c(G,data), (4.33) 

ll-llL 2 (o,f;^, 2 (n)) < c(G, data). (4.34) 

Similar estimates for d hold. With these estimates in hand, we can go back to equation 
(I4.14p for V = d t u and in a similar way to conclude 

l|ut||i,oo( 0) f ;W n,2(Q ii? 3)) + ||utt|| i 2( ,f;i, 2 (n,.R 3 )) < c(G, data), (4.35) 

ll^llL°°(o,T;iy 3 > 2 (n,i?3)) + 1 1 iit 1 1 L2(p,f;W 2 < 2 (n,R?)) ^ c(G, data). (4.36) 
As for the final estimates 

l|ti|lL 2 (o,f ; vi/ 4 > 2 (c,R 3 )) + l|p|lL°°(o,f;iy 3 > 2 (n)) — c(G, data), (4.37) 

one needs to combine with the transport equation. The treatment is similar to Section 
14.31 and we omit the details. 

Summarizing the previous considerations, we are allowed to state the following result. 

Proposition 4.1 Let T > be given. Suppose that {§, u} is a strong solution of the 
Navier-Stokes-Fourier on the time interval [0, T], T < T, the existence of which is claimed 
in Proposition ^ . 1\ Assume, in addition, that 



sup ||V x .fi(V)||L~(fv? 3 ) < G. (4.38) 
te[o,f] 



Then 



sup 

te[o,f\ 



U(t, -)\\w^(Q) + \W(t,-)\\w^(Q) + ||u(t, Ollw/ 3 . 2 ^ 3 )] < c(G,data,T). (4.39) 
In particular, if T is the maximal existence time, then T = T. 
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5 Conclusion 



To begin, the standard parabolic regularity theory implies that the solutions belonging to 
the class specified in Proposition 13.11 are, in fact, classical solutions, meaning all relevant 
derivatives appearing in the system (11.21 - fl~4l) are continuous in the open set (0, T) x Q, 
the fields g are continuously differentiable up to the boundary dfl, and the bound- 
ary conditions ( 11.71) . ( II. 8p hold, see Matsumura and Nishida [21], [22J. Consequently, 
combining Proposition 13. II with Proposition 14. 1 1 we obtain the conclusion of Theorem 12. 11 
We conclude the paper by comments concerning the hypotheses of Theorem 12.11 In 
the light of the results obtained by Fan et al. (8], the "optimal" regularity criterion is 
expected to be the bound 



r iivxui 

Jo 



L°°(n;i? 3x3 ) 



dt < oo. 



Note, however, that, unlike [8], we have to handle the transport coefficients that depend 
effectively on the temperature. The advantage of this approach is that, again unlike [8], 
the Navier-Stokes- Fourier system is known to possess a global in time dissipative solution. 

The regularity of the strong solutions seems optimal, at least in view of the result 
[T2] . Indeed the time derivative d t u of the strong solutions should be at least of class 
L 2 (0,T; L°°(fi; R 3 )) for the technique of [12] to be applicable. 

The next remark concerns the possibility to include the more general class of Newto- 
nian stress, namely 

S(0, V x u) = (v x u + V^u - jjdiv^ul) + ^(^div^ul, 

where we have added the bulk viscosity, with the coefficient rj = ^(fl). Clearly, we can 
show the same result if 

— — = A — a constant. 
Moreover, our method still applies if 

W) 

where small is determined in terms of the optimal constant A in the elliptic estimate 



"small" . 



Av H — div^vl 
3 



> A||v||vi/2,2m ; R3), v|an = 0. 

L 2 (Q;i?3x3) 



The general case may be handled in the following way: 
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1. After establishing the uniform bounds in Section 13.11 apply the theory of Krylov 
and Safonov [IT] for parabolic equation in n on- divergence form to equation (13.11) to 
obtain uniform bounds on $ in the Holderian norm. 

2. Use the regularity theory for elliptic systems with Holder coefficients in Section I4T31 

Finally, we remark that the methods of the present paper could be extended to other 
types of boundary conditions of the Navier type for the velocity and to more general 
classes of domains Q including certain unbounded domains like exterior domains or a 
half-space. 
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